Current induced torques and interfacial spin-orbit coupling: Semiclassical Modeling 
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In bilayer nanowires consisting of a ferromagnetic layer and a non-magnetic layer with strong 
spin-orbit coupling, currents create torques on the magnetization beyond those found in simple 
ferromagnetic nanowires. The resulting magnetic dynamics appear to require torques that can be 
separated into two terms, damping-like and field-like. The damping-like torque is typically derived 
from models describing the bulk spin Hall effect and the spin transfer torque, and the field-like 
torque is typically derived from a Rashba model describing interfacial spin-orbit coupling. We 
derive a model based on the Boltzmann equation that unifies these approaches. We also consider an 
approximation to the Boltzmann equation, the drift-diffusion model, that qualitatively reproduces 
the behavior, but quantitatively fails to reproduce the results. We show that the Boltzmann equation 
with physically reasonable parameters can match the torques for any particular sample, but in some 
cases, it fails to describe the experimentally observed thickness dependences. 



I. INTRODUCTION 

Spintronic applications like spin-transfer-torque mag- 
netic random access memory (STT-MRAM) or magnetic 
domain wall-based devices require advances in materials 
to reach their full potential. The goal of improving these 
materials has led to the study of bilayers consisting of 
ferromagnetic layers and non-magnetic layers with strong 
spin orbit coupling. Recent measurements on such sys- 
tems have demonstrated efficient switching of magnetic 
tunnel junctions; 1 - like those used in STT-MRAM, and 
efficient current-driven domain wall motion^ 

There are a number of physical processes- 3 in these 
systems that contribute to the magnetization dynam- 
ics as described by the Landau-Lifshitz-Gilbert equation. 
These include the typical micromagnetic contributions, 
like interatomic exchange, magnetostatic interactions, 
magnetocrystalline anisotropy, and damping, as well as 
the adiabatic and non-adiabatic spin transfer torques 4 ^— 
that are typically added to account for the coupling be- 
tween the magnetization and the electrical current flow- 
ing through it. In the bilayers of interest here, there 
are additional contributions that have received extensive 
attention. These arise from the spin-orbit coupling in 
the non-magnetic layer and from the enhanced spin-orbit 
coupling at the interfaces between layers. 

These additional contributions have been modeled in 
terms of two different pictures. One picture 1 ^ assumes 
that the layers are thick and the two layers have their 
bulk properties. A current flowing through the non- 
magnetic layer with strong spin-orbit coupling generates 
a spin current perpendicular to the interface (the spin 
Hall effecti^ 3 - - — ) When this spin current impinges on the 
interface, there is a spin transfer torque 1 ^— on the mag- 
netization of the magnetic layer. The details of the torque 
in this picture are determined by the bulk spin Hall angle 



in the material with strong spin-orbit coupling and the 
mixing conductance. The other pictured— assumes two 
dimensional transport that can be described by a Rashba 
model, similar to those used to describe spin-orbit cou- 
pling in two-dimensional electron gases ^ The Rashba 
model gives direct coupling between the magnetization 
and the flowing current. Both models give qualitatively 
similar results, that is torques along the M x (jxz) and 
M x [M x (j x z)] directions, where M is the magneti- 
zation, j is the in-plane current density and the interface 
normal is in the z direction. We refer to the first torque 
as a field-like torque because it has the same form as pre- 
cessional torque around an effective field in the —jxz 
direction. The second torque has the same form as a 
damping torque toward a field in that same direction and 
we refer to it as a damping-like torque3i 

Both models have strengths and weaknesses. The 
Rashba model treats the strong spin-orbit coupling at 
the interfaces between the materials but treats the trans- 
port as two-dimensional. The layer thicknesses are usu- 
ally comparable to mean free paths and spin diffusion 
lengths requiring a three dimensional description of the 
transport. On the other hand, the spin-Hall-cffcct spin- 
transfer-torque model treats the three dimensional as- 
pect of the transport, but ignores any contributions from 
the modification of the spin-orbit coupling near the in- 
terface. The non-magnetic layer and the magnetic layer 
affect the electronic structure of each other close to the 
interface and the interaction can significantly change the 
spin-orbit coupling there. In particular, the proximity 
to the ferromagnet can induce a moment in the material 
with strong spin-orbit coupling giving a thin layer where 
the magnetism and the spin-orbit coupling coexist^ 

Attempts to develop predictive models face the com- 
plication that the experimental structures deviate signifi- 
cantly from the ideal structures treated theoretically. Ex- 
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perimental indication o 28 ' 29 that interfaces of Co grown on 
Pt have different properties than interfaces of Pt grown 
on Co argue strongly that the details of the interface 
structure are both nontrivial and important. Unfortu- 
nately, the interfaces are not well enough characterized 
to know what types of disorder might be present. There 
may be significant interdiffusion at the interfaces be- 
cause, for example, Pt alloys with Co in the bulk. There 
is also significant lattice mismatch between the materials. 
This mismatch could promote thickness fluctuations and 
dislocation formation. Without measurements of atomic 
scale structure of the experimental samples, it is impos- 
sible to know how important such defects are to the be- 
havior of the system. 

Motivated by the uncertainty in the details of the ex- 
perimental structures and the goal of incorporating the 
strengths of existing models, we develop simple semiclas- 
sical models for these systems. One approach is based 
on the Boltzmann equation and the other on the drift- 
diffusion equation. Both capture the essential physics of 
the models that have been used so far and provide a test 
for whether a model based on bulk properties and en- 
hanced spin-orbit coupling at the interface can account 
for the experimental behavior. We find that these mod- 
els are general enough to reproduce the torques mea- 
sured in any single sample for reasonable values of the 
parameters, but not all samples with a single parame- 
terization. In Sec. |TT] we summarize the experimental 
evidence for the different interpretations. We also sum- 
marize the micromagnetic simulations that provide the 
basis for these interpretations. We briefly present the 
semiclassical transport models in Sec. lIIII and then apply 
them in Sec. IIVI to the model for a bulk spin Hall cur- 
rent in the non-magnetic layer leading to a spin transfer 
torque at the interface with the ferromagnet. In Sec. |V] 
we add in the interfacial spin orbit coupling and show 
that this captures the important physics that is included 
in the Rashba-model approaches. 



II. EXPERIMENTAL RESULTS AND 
THEORETICAL IMPLICATIONS 

Recent experiments on multi-layer structures report 
evidence for current-induced torques due to spin-orbit 
coupling. These torques are reported to be large enough 
to modify the magnetization dynamics and may be uti- 
lized to facilitate spintronic applications. Various experi- 
ments report conflicting results for the size and direction 
of the torque. Many of these values are inferred from 
measurements of the resulting dynamics in conjunction 
with simulations. This section aims to summarize the 
evidence from experiments, in conjunction with simula- 
tions, for different forms of the torque. 

The recent growth of interest in bilayer systems began 
with a series of experiments by Liu et alr^^— and Miron 
et a/^ i 33 ! 34 Both sets of experiments treat systems of a 
substrate, non-magnetic layer, ferromagnetic layer, and 
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FIG. 1: (color online) Schematic of a bilayer structure that 
consists of a ferromagnetic layer and a non-magnetic metal 
layer with strong spin-orbit coupling. The spin Hall effect 
in the non-magnetic layer bends electron trajectories and in- 
jects electrons with proper spin direction into the adjacent 
magnetic layer, thereby generating the damping like torque. 
In this illustration, the spin Hall angle for the non-magnetic 
layer is assumed to be positive. The figure also highlights the 
interface region between the two layers, where the Rashba 
spin-orbit coupling may have sizable magnitude. 

capping layer, see Fig.[TJ although the details of each dif- 
fer from experiment to experiment. The authors of the 
first set of experiments interpret their results in terms 
of a dominant damping like torque that they attribute 
to the spin Hall effect and as subsequent spin transfer 
torque. On the other hand, the authors of the second set 
of experiments interpret their results in terms of an im- 
portant field like torque that they attribute to interfacial 
spin orbit coupling. 

In the case of the spin Hall effect, electron trajecto- 
ries are preferentially scattered in different directions de- 
pending on their spin directions. For instance, in a non- 
magnetic material with positive spin Hall angle, elec- 
trons are scattered more strongly into directions such 
that (v; x Vf) • S is positive, where vvn is the electron 
velocity before (after) the scattering and S is the elec- 
tron spin direction. Thus in the bilayer system in Fig. [T] 
the spin Hall effect in the non-magnetic layer injects elec- 
trons with spin along +y direction into a magnetic layer 
if j is along +x direction and the spin Hall angle is pos- 
itive. As it does for perpendicular flow of electrons in 
multilayer systems, this spin current causes spin transfer 
torquesji£~— the injected electrons giving rise to a spin 
torque along the direction — M x (M x m), where m is 
the direction of the magnetic moment carried by injected 
electrons and points along — y direction since m and S 
are anti-parallel. This torque amounts to the damping- 
like torque — M x [M x ( j x z)] arising from the spin Hall 
effect. 

Liu et al. examine various magnetization dynam- 
ics caused by the spin Hall effect. In Pt/Py bilayer^ 
they use the spin-torque ferromagnetic resonance tech- 
nique to show that the spin Hall effect is strong enough 
to cause magnetic precession. Through the resonance 
line shape analysis, they quantify the magnitude of 
the spin-Hall-effect-induced damping like torque and the 
torque due to the Oersted field. They estimate the spin 
Hall angle of Pt to be about +0.076, which is about 
two orders of magnitude larger than the correspond- 
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ing value in n-doped GaAs. 35 ' 36 This demonstrates that 
the spin Hall effect can be a realistic tool to enhance 
spin torque efficiency. In their subsequent work 1 - for the 
Ta/Co4oFe4oB2o bilayer, they demonstrate that the spin- 
Hall-induced damping-like torque can switch the magne- 
tization in a reliable and efficient way, facilitating the 
development of magnetic memory and nonvolatile spin 
logic technologies. From three different measurements 
(spin-torque ferromagnetic resonance, current-dependent 
anomalous Hall effect, threshold current density for mag- 
netization switching), they estimate the spin Hall angle 
of Ta to be —0.12 to —0.15. This angle is of opposite sign 
compared to Pt but most notably, larger by about a fac- 
tor of two. In their still later work, 31 they demonstrate 
that the magnetization switching can be achieved by the 
spin Hall effect in Pt(2 nm)/Co(0.6 nm)/A10a; systems 
as well, although the spin Hall angle in Pt is smaller. 
Most recently, they have determined an even larger spin 
Hall angle in Wj2£ 

On the other hand, Miron et a^ ' 33 ' 34 explore effects of 
the Rashba spin-orbit coupling with focus on the system 
Pt(3 nm)/Co(0.6 nm)/A10a;(2 nm). The layer structure 
of the film breaks inversion symmetry and gives rise to 
perpendicular magnetic anisotropy strong enough that 
the ground state for the magnetization is out of plane. In 
these experiments, the nominally symmetric control sys- 
tem Pt(3 nm)/Co(0.6 nm)/Pt(3 nm), shows much weaker 
effects, leading to the inference that the broken inversion 
symmetry for the layer plays a crucial role. 

Miron et al interpret their results in terms of a large 
field- like torque as expected from the Rashba model^ 5 . for 
interfacial spin-orbit coupling. Calculations 2 ^— predict 
an effective field 3 ^ 



when a system is subject to the Rashba spin-orbit cou- 
pling of the form aR.(k x z) • cr. Here M s is the saturation 
magnetization, P is the spin polarization, and /xb is the 
Bohr magneton. The effective field Hr generates the 
field-like torque — 7M x Hr cx — M x (z x j), where 7 is 
the gyromagnetic ratio. 

In Ref. [33, the authors measure the reversal of the 
perpendicular magnetization. They find that transverse 
in-plane magnetic fields enhance the nucleation and that 
this enhancement is modified by currents flowing through 
the film. In-plane currents j flowing through the sys- 
tem enhance (suppress) the nucleation of reversed mag- 
netic domains when the direction z x j is parallel (anti- 
parallel) to externally applied in-plane magnetic fields. 
In the nominally symmetric control sample, the effect 
of in-plane current is much weaker. The authors inter- 
pret the effect of in-plane currents as evidence for the 
current-induced effective field along z x j, as predicted 
theoretically for systems with Rashba-like spin orbit cou- 
pling. Experimentally^ the magnitude of the effective 
field is proportional to j with a proportionality constant 
(1.0 ±0.1) x 10~ 12 T/(A- m~ 2 ). This value implies a R « 



0.2 eVnm^ which is comparable to or = 0.3 eVnm 
reported for Bi/Ag(lll) surface 4 ^ and comparable to 
a recent first principles calculation result for a Pt/Co 
bilayer 4 1 - 

Further evidence for a field-like torque is found in 
measurements 2 - of current-driven domain wall (DW) mo- 
tion in the same system. As the driving current density 
goes up, the domain wall velocity wdw increases to speeds 
up to ss 400 m/s, more than three times faster than previ- 
ously measured current-driven domain wall velocities 42> 
The authors claim that this velocity is twice as large 
as the rate v s — |j|(Pg^B)/(2eMs) of the spin angu- 
lar momentum transfer, where g(~ 2) is the gyromag- 
netic ratio. Even up to such high domain wall speeds, 
the domain walls apparently did not undergo structural 
instability (Walker breakdown 1 ' ) According to con- 
ventional theories 6 -^ of current-driven domain wall mo- 
tion, domain wall velocities above v s are possible below 
the breakdown current density jw- However, as the ra- 
tio j)dw/"s increases, the breakdown current density de- 
creases such that conventional theories cannot explain 
the velocities measured in this experiment i 2 - 

As a possible explanation, the authors suggest 2 - that 
the current-induced effective field Hr could increase the 
breakdown current density. Consider the two low en- 
ergy structures [Figs.[2fa) and (b)] of the Bloch domain 
wall. These walls differ from each other because they 
have opposite magnetization directions at the domain 
wall center. Since the two directions are parallel and 
anti-parallel to Hr, the effective field either stabilizes or 
destabilizes the corresponding Bloch domain wall struc- 
tures [Fig. H^d)]. Recalling that the conventional spin 
torques produced by the current tend to shifts <fi away 
from the low energy values ((f) — and it), the deeper 
energy valley implies an enhanced threshold current den- 
sity to escape the valley. Since the Walker breakdown oc- 
curs when ip cannot remain stationary, the Walker break- 
down threshold current density becomes the larger of the 
two, which is larger than the Walker breakdown thresh- 
old value without Hr. 

Another interesting feature of the experiment 2 is that 
contrary to previous experiments ) 42 ' 45 ' 46 which report do- 
main wall motion along the electron flow direction, Ref. @ 
reports domain wall motion against the electron flow. 
Within the scope of the conventional theories* 6 ^ the re- 
versed domain wall motion implies negative P/3, where P 
is the polarization of the current and j3 is the dimension- 
less coefficient of the non-adiabatic spin transfer torque. 
All of the qualitative feature of the experiment 2 would be 
explained if there were an Hr that gives Walker break- 
down threshold enhancement, negative P/3 for reversed 
motion, and \/3/a\ > 1 for velocity enhancement (a is the 
Gilbert damping constant). 

Further measurements^ 4 on the same Pt(3 
nm)/Co(0.6 nra)/A10 x (1.6 nm) system reveal a 
problem in this simple theoretical picture. Based on 
measurements of bipolar switching in a tilted field, 
the authors conclude that the switching is due the 
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FIG. 2: (color online) Two possible structures of the Bloch 
domain wall. In (a), the magnetization at the center of the 
domain wall points along zxj and in (b), it points along —zxj. 
Panel (c) shows the domain wall energy as a function of the 
domain wall tilting angle <f> without the effective field Hr. 
Panel (d) shows the energy with the effective field assuming 
that qrP is positive and Hr prefers <j> = Q. 



additional presence of a damping-like torque along 
M x [M x (z x j)]. The damping-like torque is either 
parallel or anti-parallel to the nonadiabatic STT at the 
domain wall center and thus modifies the domain wall 
velocity] 47 ! 48 A recent theoretical study 4 ^ demonstrates 
that high «dw against electron flow direction is possible 
(Fig. [3]) even when both P and j3 are positive and /3/a is 
smaller than 1, if both the field-like torque — M x (zxj) 
and the damping-like torque M x [M x (zxj)] are 
sufficiently large. 

Miron et alS^- estimate that the spin Hall 
contribution 31 to the damping-like torque is not 
strong enough to explain the bipolar switching and con- 
clude that the damping-like torque arises mainly from 
the Rashba spin-orbit coupling, based on the observation 
that the efficiency of the bi-polar switching increases 
with the magnetic anisotropy of the cobalt layer and the 
oxidation of the aluminum layer. Calculations 3 - ' 47 ' 49 ! 50 
confirm that the Rashba spin-orbit coupling can give rise 
to the damping-like torque. Note however that the spin 
Hall contribution and the Rashba spin-orbit coupling 
contribution to the damping-like torque have exactly the 
same structure, making the distinction difficult. In this 
regard, recent micromagnetic calculations 47 show that 
the spin Hall effect and the Rashba spin-orbit coupling 
contributions give qualitatively different domain wall 
motion if there is no field-like torque as might be 
expected if the spin Hall effect is the dominant source of 
the effect. Without the field- like torque, the domain wall 
motion at high current densities is along the electron 
flow direction (if f3P > 0) whereas with the strong 
field-like torque, the domain wall motion against the 
electron flow is possible even for j3P > (solid line in 
Fig. |3l). 

In contrast to Miron et al.'s interpretation^ Liu 
et at 31 conclude that the damping-like torque in 
Pt(2.0 nm)/Co(0.6 nm)/Al(1.6 nm) arises mainly from 
the spin Hall effect in Pt and the Rashba spin-orbit cou- 
pling contribution is negligible. As evidence for the lat- 
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FIG. 3: (color online) Terminal domain wall velocity i>dw 
as a function of the current density j = jx. The dotted 
line gives the prediction of the Landau-Lifshitz-Gilbert equa- 
tion when the current generates only the conventional (adi- 
abatic and nonadiabatic) spin transfer torques and the solid 
line the prediction when the current also generates field-like 
and damping-like torques. For this calculation, we consider 
a nanostrip of length x width x thickness = 2000 nm x 
200 nm x 5 nm, with perpendicular magnetic anisotropy and 
the materials parameters as follows: 7/(2-71-) = 28 GHz/T, 
M s = 1 x 10 6 A/m, A cx = 1.3 x 10" 11 J/m, P = 0.7, 
K u = 1.5 x 10 6 J/m 3 , a = 0.5, /3 = 0.25, and a R = 
0.7 x 10~ 10 eVm. For the conventional calculation (dotted 
curve), the non-linearities are due to Walker breakdown. In 
the solid curve, the small blip at a small negative current is 
caused by the chirality switching of domain wall due to the 
spin-orbit-related field-like torque. 



ter conclusion, they report that the ratio between the 
current-induced effective field along zxj and the in-plane 
current density is more than 75 times smaller than the 
corresponding ratio in Ref. l33l . A third experiment on 
a similar system Pt(3.0 nm)/Co(0.6 nm)/A10 a; (1.8 nm) 
found that this ratio is about 3.4 times smaller than the 
ratio reported in Ref. |H. 

Experiments 5 ^ that measure the field-like torque on 
Ta(1.0 nm)/Co4oFe 4 oB2o(1.0 nm)/MgO(2.0 nm) find 
that the ratio of the effective field to the current is about 
23 % of the corresponding ratio in Pt/Co/AlOa; as re- 
ported in Ref. l33l . Interestingly, the current-induced ef- 
fective field in the Ta/CoFeB/MgO system is along —zxj 
and thus opposite to the corresponding direction in the 
Pt/Co/AlOa; system if the Ta and MgO layers in the 
Ta/CoFeB/MgO system are matched with Pt and AlOa; 
layers in the Pf/Co/AlO^ system, respectively. It was 
also reported that the torque to current ratio decreases 
by more than an order of magnitude when the thickness 
of the Co4oFe4oB 2 o layer in the Ta/CoFeB/MgO system 
increases slightly from TO nm to 1.2 nm. 

The thickness dependence was investigated 
systematically^ 3 - for the two wedge systems, 
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Ta((iTa)/Co2oFe6oB2o(l nm)/MgO(2 nm) and 
Ta(l nm)/Co 20 Fe 6 oB 20 (tcoFeB)/M g O(2 nm). When 
the Ta layer thickness c?Ta changes by 1 nm, the effective 
field along the z x j direction changes its magnitude by 
nearly two orders of magnitude. For small rfxa 0.6 nm, 
the sign of the effective field is opposite to that in 
the larger dxa regime and agrees with the sign of the 
Pt/Co/AlOj; system^ The damping-like effective field 
along the (z x j) x M direction is also sensitive to dx a 
with the sign change at dxa ~ 0.5 nm. The sign change 
of the damping-like effective fields is interpreted as 
an evidence of competition between Rashba spin-orbit 
coupling and spin Hall effect. 1 

Measurements of current-driven domain wall motion 
in multi-layer structures containing non-magnetic heavy 
metal layers and ferromagnetic layers but without ox- 
ide layer o 29 ' 54 are reported. The motion is also very 
sensitive to layer thicknesses.— The domain wall veloc- 
ity can be up to almost 1 km/s in certain multi-layer 
structures with the Pt layer thickness of about 1 nm but 
changes very quickly as the Pt thickness varies. Inter- 
estingly, very high domain wall speed (~ 1 km/s) are 
observed only when the domain wall moves against the 
electron flow, implying that the origin of the reversed do- 
main wall motion is probably correlated with the mecha- 
nism behind very high domain wall speed. Several other 
experiment s 28 ' 55 " — also report reversed domain wall mo- 
tion in ultrathin multi-layer systems containing Pt layers. 

The analysis of domain wall motion described above 
considers four current induced torques, the the adiabatic 
and non-adiabatic spin-transfer torques and two torques, 
damping-like and field-like, that depend on the layer 
structure of the device. The first two torques depend on 
the gradient of the magnetization but are determined by 
bulk properties. The last two are independent of the gra- 
dient of the magnetization. Other possibilities that de- 
pend on the gradient of the magnetization and the layer 
structure are allowed by symmetry^ and may be impor- 
tant for the dynamics. In addition, recent calculations 58 
suggest that a current-independent torque due to the 
Dzyaloshinskii-Moriya interaction 59 might provide an al- 
ternate mechanism for stabilizing a moving domain wall 
above the nominal Walker-breakdown field. In this pa- 
per we only consider the current-induced torques that are 
independent of the gradient of the magnetization. 



III. SEMICLASSICAL MODELS 

To explore possible mechanisms for the torques oper- 
ative in these systems, we develop semiclassical mod- 
els that allow for easy exploration of parameter space. 
We use a Boltzmann equation approach and the sim- 
pler drift-diffusion approach. The Boltzmann equation 
is better suited to describe in-plane transport but the 
drift-diffusion approach is simpler and provides a useful 
language to describe the physics. 

The Boltzmann equation approach developed by Cam- 



ley and Barnas^S is the simplest model that describes 
current-in-the-plane (CIP) giant-magnctoresistance 
(GMR). In the Boltzmann equation, the variables are 
the distribution functions. The distribution function 
accounts for electrons moving in all directions even 
though the total current only points in a single direction. 
This generality allows the approach to describe the flow 
of spins between the layers even though the current flows 
in the plane of the interfaces. 

The drift-diffusion approach of Valet and Fert 61 
is based on integrating the distribution function in 
the Boltzmann equation to derive transport equations 
that depend on the densities and currents. It has 
had wide success describing current-perpendicular-to- 
the-plane (CPP) GMR, but does not describe CIP GMR. 
It fails because it does not describe the flow of spin cur- 
rents between the layers when the current flows in the 
plane of the layers. However, this limitation may be less 
important in the bilayer systems of interest here. In ma- 
terials with strong spin-orbit coupling, like Pt, spin cur- 
rents do flow perpendicular to the charge current because 
of the spin Hall effect^ so that the drift-diffusion ap- 
proach does qualitatively describe the physics. However, 
we compare the two approaches below, and show that 
the drift-diffusion approach differs quantitatively from 
the Boltzmann equation for similar reasons to its quali- 
tative failure for CIP GMR. 

Reference [63| describes the matrix Boltzmann equation 
we use in this paper. It is a generalization of the model 
used to describe CIP GMR^S and is based on a very 
simplified model for the electronic structure. We treat 
all Fermi surfaces as spherical and as having the same 
Fermi wavevector. This approach ignores the details of 
the Fermi surfaces, which are undoubtedly important for 
specific systems, particularly when including spin-orbit 
coupling. However, the scattering mechanisms are both 
unknown and uncharacterized, so we feel that it is ap- 
propriate to consider models in which scattering rates 
and other physical processes are parameterized and the 
details of the electronic structure are neglected for sim- 
plicity. By performing appropriate integrals over the dis- 
tribution function, the Boltzmann equation can be trans- 
formed into a drift-diffusion equation like that given in 
Ref. 164 . The parameterized processes in the Boltzmann 
equation then have a simple connection to those in the 
drift-diffusion equation. 

One such process that we include through parameter- 
ization is the spin-dependent conductivity in the ferro- 
magnetic layer. We model the spin dependence by using 
spin-dependent scattering rates. Such spin-dependent 
scattering is physically sensible as it is believed to play a 
bigger role in the polarization of the current in these ma- 
terials than the electronic structure itself. For example, 
in Co and Ni, the current is expected to be dominated 
by majority carriers, even though there are more minor- 
ity carriers at the Fermi surface. In the drift diffusion 
model, the spin-dependent scattering leads to a differ- 
ent conductivity for the majority electrons than the 
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minority electrons cr*. This difference is parameterized 
in terms of the spin polarization of the current, defined 
through P = (o 4 " - (T l )/{^ + o+). The drift-diffusion 
transport equations in the ferromagnet are 

j = aVu- P<tV(M- /x s ), 



continuity equations in the ferromagnet are 



2e 



2e 



(2) 
(3) 



where j is the charge current density, and Q is the tensor 
spin current density where the first index is the spatial 
component and the second index the spin component, 
/i is the electrochemical potential, such that negatively 
charged electrons diffuse against the gradient giving an 
overall positive sign for the first term in Eq. ([2]). Similarly 
fi s is the spin chemical potential, which is a vector along 
the direction of the spin accumulation, and the unit vec- 
tor M is the direction of the magnetization. The minus 
sign in the second term of Eq. (|2|) arises because majority 
electron spins are aligned opposite to the magnetization. 
These two signs are typical of possible sources of confu- 
sion in this subject matter. They arise because with the 
charge on the electron is negative and angular momenta 
and moments are in opposite directions. 

The equality of the Fermi surfaces would also allow 
for perfect transmission of electrons across the interface 
between the materials. In the Boltzmann equation, we 
include spin-dependent reflection by the addition of a 
spin-dependent sheet potential (delta function) at the 
interface. Choosing the strength of this delta function 
allows us to tune the spin-dependent interface resistance 
to any arbitrary valued In the drift-diffusion approach, 
the spin-dependent reflection becomes a spin-dependent 
interface resistance or conductance as used in the closely 
related circuit theory^ 

While the overall structure of the Boltzmann equation 
approach is the same as that published earlier, there are 
some modifications. One difference is the treatment of 
dephasing. In a ferromagnet, spins on different parts of 
the Fermi surface precess at different rates and travel 
with different velocities. These differences, combined 
with scattering between different parts of the Fermi sur- 
face, cause the processing spins to rapidly become out of 
phase with each other^s In Ref. |63|, the transverse spin 
accumulation and current are forced to zero in the fer- 
romagnet to account for this dephasing of the transverse 
spin population. Here, we allow for transverse spin accu- 
mulation in the ferromagnet but build in rapid spin pre- 
cession and explicitly account for the processes that cause 
dephasing. The simplified model of the Fermi surfaces 
that we use can lead to underestimation of dephasing 
processes. We have tested this approximation by adding 
an explicit dephasing term. While such a term quantita- 
tively changes the spin accumulation in the ferromagnet, 
we find that it does not change the calculated torques. 

The same dephasing process is absent in a drift- 
diffusion model but can be included by adding an explicit 
dephasing term. In this approximation, the steady-state 



Vj = 



1 



1 



ViQij = (s x M)j 



' ; v r sf 
M x (s x M) 



1 

T~dp 



(4) 



(5) 



where the spin accumulation s is proportional to the spin 
chemical potential s = Af s fi s and the precession time 
Tex = h/A is related to the exchange splitting between 
the magnetization and the spin accumulation. Repeated 
indices are summed over. The first term on the right 
hand side of Eq. ([5]) is the precession in the exchange 
field, the second term is the spin flip scattering that re- 
duces all components of the spin accumulation, and the 
last term is the dephasing that reduces only the parts of 
the spin accumulation transverse to the magnetization. 
Setting the transverse spin accumulation to zero, as done 
in earlier Boltzmann equation calculations 53 and in mag- 
netoelectronic circuit theory^ is equivalent to taking the 
limit that the dephasing time goes to zero. 

Both the size of the spin Hall effect and its underly- 
ing mechanism are controversial. The theory for the spin 
Hall effect is related to that for the anomalous Hall effect, 
a subject that has been controversial for decades^ Mea- 
surements of the spin Hall angle (the ratio of spin Hall 
and charge conductivities) for various materials span a 
range of values. Part of the variation may result from the 
sensitivity of the extraction of the spin Hall angle from 
experimental data to other material parameters needed 
to model the experiments^ Measurements^ show that 
the spin Hall effect in Ta is bigger and of the opposite sign 
of that in Pt, in agreement with previous calculations. 70 
The agreement between these trends in theory and exper- 
iment argues for an intrinsic origin of the effect. How- 
ever, the calculated spin Hall conductivity for Pt appears 
to be approximately an order of magnitude too small in 
comparison to the measured value. Spin Hall angles of 
approximately the right order of magnitude have been 
computed 71 for the extrinsic contributions of various im- 
purities in Cu and Au. 

With this uncertainty in the mechanism for the spin 
Hall effect, we use the form of scattering appropriate for 
the extrinsic skew scattering contribution for computa- 
tional simplicity. In the Boltzmann equation, we include 
skew scattering as described by Engel et ali2& but gen- 
eralize their results to include scattering that leads to 
the inverse spin Hall effect in addition to the scattering 
that gives rise to the spin Hall effect. These scattering 
terms connect the current with a perpendicular spin cur- 
rent and vice versa. Both our approach and the earlier 
work^li neglect the scattering processes that couple spin 
currents to spin currents moving in other directions. Such 
process contribute to spin relaxation, which we include 
as a phenomenological spin flip scattering process. De- 
tails of the scattering and the default parameters we use 
are given in Appendix 
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In the non-magnetic material, the explicit forms of the 
charge and spin currents in the drift-diffusion approxi- 
mation we use are 6 ^ 



j = <jVh - ct sh (V x fj, s ), 



2e 



2c 



O-SHCyfeVfe^. 



(6) 
(7) 



where a is the conductivity, ctsh is the spin Hall con- 
ductivity coupling the spin and charge currents to the 
charge and spin potentials, and eyfe is the Levi-Civita 
symbol. As with the Boltzmann equation, we neglect a 
term corresponding to the spin Hall effect coupling the 
spin current to the spin potential, assuming that it is 
small. 

In both models, the torque on the magnetization is 
given by the torque between the magnetization and the 
spin accumulation 
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-M x s 



7 



TdpM 



M x (M x s), 



(8) 



' where the gyromagnetic ratio 7 = g^/Ti converts from 
angular momentum (spin density) to magnetization (so 
T is a term in the Landau-Lifshitz-Gilbert equation^). 
The second term, which is not present in the Boltzmann 
equation calculations, captures the torque due to the de- 
phasing of the electron spins as they precess in the ex- 
change field. 

If there is no coupling of angular momentum into the 
lattice (spin-orbit coupling or spin-flip scattering) it is 
straightforward to relate this torque, Eq. ©, to the di- 
vergence of the spin current. In the ferromagnet, where 
there is spin-flip scattering but no other spin-orbit cou- 
pling, the corrections to torque being simply the diver- 
gence of the spin current can be found from Eq. ([5]) . For 
the components perpendicular (_L) to the magnetization, 
the torque is 



T = - 7 (1 - 0) (V ■ Q) ± + 7/3M x (V • Q), 



1 



r sf 1 + e 



Tdp 



7*cx 
Tat ' 



(9) 



where (V ■ Q) ± = — Mx [Mx (V-Q)] is the component of 
the divergence of the spin current that is perpendicular to 
the magnetization. Equation ([§]) indicates that the spin 
torque is in general not simply given by the divergence of 
the spin current but possesses an additional component 
Mx VQ arising from the presence of spin relaxation. 
For the parameter set used here, (3 is negligible. 

An important difference with the previously 
published 63 formalism for the Boltzmann equation 
is the inclusion of spin-orbit coupling at the interface. 
This is done by including an additional term in the 
interface potential 



V(r) 



h 2 k F 



S(z) uq + u cx cr ■ m + mrct • (k x z) , (10) 



where the interface is in the z direction at z = 0, uq 
is the spin-independent part of the potential, u cx is the 
spin-dependent part of the potential that gives rise to 
spin-dependent reflection, Mr is the Rashba contribution, 
with k being the wave vector of an electron scattering 
from the interface, k F is the Fermi wave vector, and m 
is the electron mass. This additional term captures the 
form of spin-orbit coupling that is allowed for the sim- 
ple electronic structure assumed here. For more realistic 
band structures, the form would be much more compli- 
cated. Unfortunately, it is difficult to compare Mr with 
the qr, used in previous publications. Doing so requires 
a procedure for reducing the Hamiltonian for a three di- 
mensional system to one for a two-dimensional system. 

In Eq. (fT0|) , the last two terms can be combined to 
give a wave vector dependent field direction u(k) and 
strength u c ff(k) such that u c ff(k)u(k) = u cx rh + urIc x z. 
With respect to this direction, the majority and minority 
transmission and reflection amplitudes are 



T = 
R = 



ik z /k F - {u ± u c g) 

Up ± Meg 

ik z /k F - (u ± u cS ) 



(11) 
(12) 



Since both the magnitude and phase of the transmission 
and reflection amplitudes are different for the majority 
and minority spin components, an electron spin oriented 
along some arbitrary direction undergoes a finite rotation 
when transmitted or reflected. A part of the torque on 
the electron spin is due to the coupling to the exchange 
field and a part due to the spin-orbit coupling (Rashba 
contribution). The reaction torque on the magnetization 
can be computed from the exchange coupling between 
the spin density at the interface and the exchange field 



7 / hk F u C x , 
S ^TF I s x M 



(13) 



where the spin density s is calculated from the incoming 
wave function and the transmission amplitudes. Since 
the potential is proportional to a delta function, the 
torque density diverges but is finite when integrated over 
a finite thickness. 

The treatment we use for the interfacial spin-orbit cou- 
pling in the Boltzmann equation does not generalize eas- 
ily to the drift-diffusion equation because there are not 
any wave vector dependent quantities in that model. It 
may be possible to define a generalization of the conduc- 
tance matrix used in the magnetoelectronic circuit the- 
ory. In typical usage, the longitudinal spin components 
couple to each other and the transverse spin components 
couple to each other, but the longitudinal and transverse 
spin components do not couple. With the Rashba inter- 
action included, all spin components would couple. 

The Boltzmann equation based approach differs quite 
significantly from the approach used in which the system 
is modeled with a two dimensional Rashba model. In our 
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Boltzmann equation approach, electron spins get kicked 
when they pass through the interface, but they spend no 
time "in" the interface. In the Rashba model, the entire 
system is the interface so the electrons (and spins) are 
"in" the interface at all times. In this case, there is a 
spin accumulation that builds up in the interface. This 
spin accumulation gives rise to the strong field-like torque 
found in these models. In spite of this difference, we 
find that the both approaches give qualitatively similar 
torques. 

The remaining difference with the previously published 
formalism for the Boltzmann equation is that the bound- 
ary conditions are different. The previous version treated 
perpendicular transport and here we treat in-plane trans- 
port. Here, the outer boundaries are perfectly reflecting, 
either diffusely, specularly, or somewhere in between. 



IV. SPIN HALL EFFECT PLUS SPIN 
TRANSFER TORQUE 

In this Section, we describe the behavior of the model 
in the absence of spin-orbit coupling at the interface. In 
this limit, the spin Hall effect in the non-magnetic layer 
generates a spin Hall current that propagates perpendic- 
ular to the interface with spins pointed perpendicular to 
both the interface normal and the direction of the cur- 
rent. When this spin current hits the interface with the 
ferromagnet, angular momentum is transfered from the 
flowing spins to the magnetization as is typical for spin 
transfer torques in magnetic multilayersj 17 ' 18 i 20 

This process is shown in Fig. [I] based on calculations 
done with the Boltzmann equation described in Sec Mil 
Parameter choices are given in Table Q] in Appendix [3] 
Fig. U shows the currents, spin currents, and spin den- 
sities for a 4 nm ferromagnetic layer coupled to a 6 nm 
non- magnetic layer with the interface at z = 0. Panel 
(a) shows the distribution through the thickness of the 
films of the current flowing in the plane of the film (in 
the x-direction) . The current is greater in the ferromag- 
netic layer because it has a higher conductivity than the 
non-magnetic layer for this choice of parameters. The 
current is suppressed close to the outer boundaries be- 
cause we assume that the scattering from those interfaces 
is completely diffuse. In fact, the ferromagnetic layer is 
not thick compared to the mean free paths, so the cur- 
rent is suppressed through the thickness of the film. The 
spin current with spins aligned with the magnetization 
(^-direction) and moving in the plane also reduced from 
the bulk value, in fact more so than the current, so the 
polarization of the current is reduced from the bulk value. 
At the interface between the two materials, the current 
is enhanced in the lower conductivity layer due to elec- 
trons entering from the higher conductivity layer, and 
the current is reduced in the higher conductivity layer. 

This modification of the current near the interface is 
not captured by a drift-diffusion model. It is one source 
of the quantitative disagreement between the models as 




5 -0.1 0.0 -0.2 0.0 
Scaled currents and spin densities 

FIG. 4: (color online) Currents, spin currents, and spin ac- 
cumulations. The left panels (a, d, and g) show the current 
density (heavy lines) j x , which is flowing in the plane of the 
sample and the spin current Q xz (lighter lines), flowing in 
the ^-direction with spins aligned with the magnetization in 
the z-direction. The dotted lines indicate the bulk values. 
All currents and spin currents are dimensionless; currents are 
scaled by the bulk current in the non-magnet and spin cur- 
rents are scaled by the bulk current in the non-magnet and an 
additional factor of h/2e. The spin densities are scaled by the 
same factor two factors and vy- The middle panels (b, e, and 
h) show the spin currents, Q zy (heavy lines) and Q zx (lighter 
lines), flowing perpendicular to the layers (z-direction) with 
spins pointing perpendicular to the magnetization, i.e. the 
x- and ^/-directions. The right panels (c, f, and i) show the 
accumulation of spin perpendicular to the magnetization, s y 
(heavy lines) and s x (lighter lines). The top panels (a, b, 
c) are for the case in which there is no interfacial spin orbit 
coupling, the bottom panels (g, h, i) for the case with inter- 
facial spin-orbit coupling iir = 0.04 and no spin Hall effect 
in the non-magnet, and the middle panels (d, e, f) for the 
case when both are present. In panels (h) and (i), the spin 
accumulations have been scaled by the indicated factors. 

discussed below. 

Panel (b) of Fig. 2] shows the two components of the 
spin current with spins aligned perpendicular to the mag- 
netization and moving perpendicular to the plane of the 
film. In the non-magnetic layer, this is due to the spin 
Hall effect. The spin current is zero at the lower bound- 
ary, which is both impenetrable and has no spin-flip scat- 
tering. It increases to close to its bulk value at the inter- 
face between the non-magnet and the ferromagnet. In- 
side the non-magnetic layer, the spin current is a compe- 
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FIG. 5: (color online) Average current density in the non- 
magnetic layer. For three thicknesses of the ferromagnetic 
layer, the average current density in the non-magnetic layer 
is shown scaled by the bulk value as a function of the thickness 
of the non-magnetic layer. 



tition between the spin Hall current and a diffusive spin 
current from the spin accumulation, seen in Panel (c), 
that builds up due to the impenetrability of the outer 
interface. Panels (d-i) show calculations with interfacial 
spin-orbit coupling included and are discussed in Sec. [V] 
At and near the interface, this spin current is converted 
into a spin transfer torque on the ferromagnet. We write 
the interfacial torque in the form 



S(z) 



2e 



r d M x (M x y) + T f M x y , (14) 



where 8(z) localizes the torque to the interface at z = 
0. The dimensionless coefficients and Tf characterize 
the "damping-like" and "field-like" contributions respec- 
tively. Other terms are possible, as we discuss in Ref. H3, 
but in the three-dimensional transport calculations we 
find these other terms to be negligible for the parameters 
we consider. 

The prefactor in Eq. (fT4"|) is based on jo, which is the 
"bulk" current density in the non-magnetic layer, that 
is jo = crisiE where E is the applied electric field. This 
choice seems to be that typically made in analyses of 
experiments even though the total current is all that 
is directly measurable. The rest of the factors convert 
from current density to magnetization torque density, 
M. This choice makes sense in analyzing experiments in 
terms of the spin Hall effect because the torque is driven 
by the current density in the non-magnet. However, for 
thin films, there can be important corrections due to the 
outer boundaries and the interface with the ferromag- 
net. These corrections are shown in Fig. [3] for a variety 
of thicknesses for the two layers. The average current 
density is reduced by the diffuse scattering assumed at 
the outer boundary of the layer, but is increased by the 
(assumed) higher conductivity of the ferromagnetic layer 
when that layer is thick enough. 

For this model, with no interfacial spin-orbit coupling, 
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FIG. 6: (color online) Torques as a function of nonmagnetic 
layer thickness. The solid curves are the full Boltzmann equa- 
tion calculation and the dashed curves give the analytic ap- 
proximation based on the drift-diffusion model and the cir- 
cuit theory (Appendix[B]). The more negative curves show the 
field-like torque Tf and those closer to zero show the damping- 
like torque t<j. For both torques, the Boltzmann results have 
been calculated for three different mean free paths (labelled 
on the damping-like torques) in the ferromagnet. 



the spin transfer torque is determined solely by the trans- 
verse spin current,— just outside the magnetic layer. 
Since neither the majority transmission probability is 
zero nor the minority reflection probability is one, some 
of the transverse spin current is reflected. The reflected 
spin current is seen in the reduction of the transverse 
spin current close to the interface. Some of the transverse 
spin current is absorbed right at the interface and some 
is transmitted into the ferromagnet. In the ferromagnet, 
spin components transverse to the magnetization rapidly 
precess as they traverse the layer, as seen in the oscilla- 
tions in panel (b). Further, different parts of the Fermi 
surface precess at different rates so they dephase as the 
traverse the layer as can be seen by the decay of the 
transverse spin current in the ferromagnet in panel (b). 

The dominant spin transfer torque arises from the ab- 
sorption of the incident transverse spin current either at 
the interface or in the ferromagnet. However, not all 
of the current is absorbed, and some is rotated into the 
x-component of the spin current on reflection. The non- 
zero reflection reduces the damping like torque and the 
rotation gives rise to a small field-like torque. These 
torques are shown in Fig. [6] as a function of the thick- 
ness of the non-magnetic layer. 

Figure |6] shows the damping-like and field- like torques 
calculated with both the Boltzmann equation approach 65 
and the drift-diffusion approach. The drift-diffusion ap- 
proach gives an analytic result, Eq. (|B2j) in Appendix 151 
That result is based on the drift-diffusion model de- 
scribed in Sec. lIIII and magnetoelectronic circuit theory^ 
for transport across the interface. Both approaches give 
the same behavior as a function of the thickness of the 
non-magnetic layer. Because the spin Hall current in the 
non-magnetic layer is suppressed when the layer is thin, 
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as seen in Fig. [4j the torque is reduced when the layer 
thickness is less than a few spin diffusion lengths, which 
for this set of parameters is £^ = 2.5 nm. 

For thick layers, the value saturates, but does not sat- 
urate to the spin Hall angle 6>sh as might be expected. 
Eq. (|B2[) shows that for the drift-diffusion model, the sat- 
uration value depends on the ratio, ^Re[g^]/aj<s. When 
this ratio is small, the saturation value is reduced from 
#SH and the Boltzmann calculation and the drift-diffusion 
calculation saturate to different values. When that ra- 
tio is large, the drift-diffusion and Boltzmann equation 
results agree. However, a large value of this ratio is not 
physically realistic for systems with strong spin-orbit cou- 
pling. The mixing conductance depends mainly on the 
area of the Fermi surface in the non-magnet (as a re- 
minder our calculations assume the same Fermi surface 
for all materials), but does so in the same way that the 
conductivity does (see the expression in Table HI so it 
is difficult to increase the ratio by changing the mixing 
conductance. It is possible to decrease the conductivity 
by increasing the non-spin-flip scattering, but this also 
decreases the spin-diffusion length. For the default pa- 
rameters we consider, see Table U the value of this ratio 
is about 0.3. 

Eq. (|B2[) also shows that the torque calculated with 
the drift-diffusion approach is independent of the details 
of the ferromagnetic layer, depending only on the mixing 
conductance. The results for the Boltzmann equation, 
for which we do not have analytic results, do depend on 
the details of the ferromagnetic layer as seen in Fig. [5] for 
different values of the mean free path in the ferromagnet. 
When the mean free path is long so that the conductiv- 
ity in the ferromagnetic layer is much greater than that 
in the non-magnetic layer, the current near the inter- 
face in the non-magnet is increased (see Fig. 0} giving 
a greater spin Hall current. Another difference between 
the approaches is that the only length scale for varia- 
tion in the drift-diffusion approach is the spin-diffusion 
length. There are many more length scales in the Boltz- 
mann equation approach (see Table and these turn out 
to play a non-negligible role when £^~Rje[g^]/aj<f is not 
large. The deviation between the results of the Boltz- 
mann equation calculations and those found from the 
drift-diffusion equation should provide a note of caution 
for the extraction of physical parameters, like the spin 
Hall angle, from comparisons between experiment and 
the drift-diffusion equation. 

We conclude that without additional spin-orbit cou- 
pling at the interface between the two materials, three- 
dimensional transport models predict a torque that is 
predominantly damping-like but has a minor field-like 
contribution. In the parameter range we have studied, 
the torque is always well described by the combination 
of these two forms. The drift-diffusion approach qualita- 
tively captures the physics but can quantitatively fail in 
physically relevant parameter regimes. 



V. INCLUSION OF INTERFACIAL RASHBA 
COUPLING 

The results in Ref. [2?] show that there is an interfacial 
region with significant spin orbit coupling and exchange 
splitting. In this section, we model that overlap region 
by adding a Rashba term to the energy at the interface, 
see Eq. (jTUJ) . We find that this additional term primarily 
leads to a field-like torque and that as long as it is not 
too strong it does not significantly modify the torques 
due to the spin Hall effect. 

Previously, this region has been treated by two- 
dimensional calculations in which the electronic struc- 
ture is modified by the competition between the Rashba 
interaction and the exchange interaction^ - — Typically, 
the Rashba interaction and the exchange interaction are 
taken to be very different in magnitude so that the Fermi 
surfaces remain essentially circular. However, the spin 
eigendirections on the Fermi surfaces are modified so that 
the non-equilibrium occupation due to a current flow give 
rise to a net spin accumulation that is not aligned with 
the magnetization. This net transverse spin density gen- 
erates a field-like exchange torque on the magnetization. 

In the Boltzmann equation approach that we use, the 
Rashba interaction modifies the boundary conditions for 
the distribution functions at the interface. The result- 
ing torque is very similar to what is found in the two- 
dimensional calculations. Depending on the details of the 
parameters, the transmission probability is, on average, 
either greater or lesser for spins aligned with z x j than for 
those in other directions. The spin density at the inter- 
face is determined by the transmission probabilities and 
the incident fluxes. The bias in the transmission proba- 
bilities favors a net spin polarization aligned with z x j, 
very similar to the behavior found in the two-dimensional 
treatments. Then, through Eq. (fT3|) . there is a field-like 
torque. 

The effect of introducing the Rashba-term on the cur- 
rent distribution is shown in Fig. 01 Unfortunately, the 
spin densities at the interface are obscured by the approx- 
imations of the Boltzmann equation. In this approach, 
we assume that electrons on different parts of the Fermi 
surface are incoherent with each other. However, the 
matching conditions for the distribution functions across 
the interface are found through coherent scattering cal- 
culations. Once the scattering states are used to con- 
struct the matching conditions, the coherence between 
the incoming and outgoing states is neglected. As a re- 
sult, for each electron the spin density at the interface is 
equal to the incident amplitude times either the trans- 
mission probability |T| 2 or |l + i?| 2 , which are the same 
because the wave function is continuous across the in- 
terface. However, immediately outside the interface, the 
incident and reflected states are no longer treated as co- 
herent, on one side of the interface the spin density is 
proportional to |T| 2 and on the other l + |i?| 2 . Since there 
are electrons incident from both sides, the spin density 
at the interface is not equal to the spin density on either 
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FIG. 7: (color online) Scaled torques as a function of the in- 
terfacial Rashba interaction. The solid curves are calculated 
with a bulk spin Hall effect and the dotted curves without. 
The non-magnetic layer is 6 nm thick and the ferromagnetic 
layer is 4 nm. The damping-like and field-like torques are 
labeled with Td and T{ respectively. The jitter is due to nu- 
merical instabilities. 



side. 

In Fig. 21 we have used a strength of the Rashba in- 
teraction that gives a torque that is comparable to that 
found from the spin Hall effect. However, this interaction 
only has a small effect on the spin currents and accumu- 
lations. We find that the spin currents for systems with 
both the spin Hall effect and the interfacial Rashba in- 
teraction can be simply and accurately approximated as 
follows. The currents in the plane, panels a, d, and g in 
Fig. [4] are essentially the same in all three systems. For 
system with both the spin Hall effect and the Rashba 
interaction, the transverse spin currents perpendicular 
to the plane, panel e in Fig. [H are a sum of the trans- 
verse spin currents found in systems with one effect or 
the other, panels b and h. The same holds true for the 
resulting torques, the torques for systems with both ef- 
fects are approximately the sum of the torques found in 
the systems with one effect or the other. 

The approximate independence of the torques due to 
the spin Hall effect and the interfacial Rashba interaction 
is illustrated in Fig. [7] as a function of the strength of the 
interfacial Rashba interaction. There are contributions 
to both the field-like and damping-like torques that in- 
crease linearly with the Rashba interaction strength up 
to large values of the torques. Eventually, the transmis- 
sion probabilities get so low, see Eq. (|T1Ij) . that all elec- 
trons are reflected and the torques go to zero. Comparing 
calculations done with and without the spin Hall effect 
shows that the interfacial coupling has very little effect 
on the torque from the spin Hall effect, particularly for 
small values of the Rashba interaction strength. With- 
out a priori knowledge of the parameters of the system, 
particularly the spin diffusion length, the spin Hall angle, 
and the interfacial spin-orbit coupling, a wide variety of 
combinations of damping-like and field-like torques are 
possible. 
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FIG. 8: (color online) Dimensionless torque components as 
a function of the thickness of the non-magnetic layer. Panel 
(a) shows the torques in the absence of the Rashba contribu- 
tion from the interfacial spin-orbit coupling, panel (c) shows 
the torques in the absence of the spin Hall effect in the non- 
magnet, and panel (b) shows the torques with both present. 
In each panel, the solid lines show the damping-like torque 
and the dotted lines the field-like torque. 



Figure [8] shows the torques as a function of the thick- 
ness of the non-magnetic layer. The torques due to the 
spin Hall effect largely depend exponentially on the non- 
magnetic layer thickness divided by the spin diffusion 
length. As can be seen from the analytic solution pre- 
sented in Appendix [Bl for very small thicknesses there 
are corrections such that the torque varies as t^ M rather 
than linearly as it would if the behavior were strictly ex- 
ponential. The exponential variation has been used^ to 
extract the spin diffusion length for particular systems. 
It is interesting to note that for this simple model, the 
sign of the product r^Tf is opposite in the two limiting 
cases, the spin Hall effect only case [Fig. EJa)] and the 
Rashba effect only case [Fig. [He)]. 

Figure ® shows the torque as a function of the ferro- 
magnetic thickness. As expected from the analytic solu- 
tion in Appendix iBl that is displayed in Fig[6j the torque 
due to the spin Hall effect only depends weakly on the 
thickness of the ferromagnetic layer. The variation is 
largely due to the variation in the current in the non- 
magnetic layer due to the presence (and variation) of the 
ferromagnetic layer. The torque due to the Rashba effect 
depends more strongly on the ferromagnetic thickness, 
but not nearly so strongly as is seen in experiment where 
changing the thickness by a single atomic layer can have 
a profound effect on the torque. 

This discrepancy suggests that the thickness depen- 
dence seen in experiment is likely due to physics beyond 
the scope of the model presented here. One possibility is 
that the strength of the Rashba interaction could depend 
sensitively on the thickness of the ferromagnetic layer. 
The sensitivity could arise from changes in the electronic 
structure of the interface or even changes in the structure 
of the interface. Since the lattice mismatch is so large, 



12 



0.03 
0.0 

-0.05 
0.03 
0.0 

-0.05 

0.03 
0.0 



(a) Spin Hall effect 



= 6 nm 



Tf 



(b) Both 



Q5 -(c) Rashba effect 



10 
/ FM (nm) 



15 



20 



FIG. 9: (color online) Dimensionless torque components as 
a function of the thickness of the ferromagnetic layer. Panel 
(a) shows the torques in the absence of the Rashba contribu- 
tion from the interfacial spin-orbit coupling, panel (c) shows 
the torques in the absence of the spin Hall effect in the non- 
magnet, and panel (b) shows the torques with both present. 
In each panel, the solid lines show the damping-like torque 
and the dotted lines the field-like torque. 



it is conceivable that the structure evolves rapidly as the 
layer is made thicker. 



VI. SUMMARY 

In this article, we have developed semiclassical mod- 
els for electron and spin transport in bilayer nanowircs 
with a ferromagnetic layer and a non-magnetic layer with 
strong spin-orbit coupling. We use a Boltzmann equation 
approach, based on a simplified electronic structure and 
also a simpler drift-diffusion model. The drift-diffusion 
framework qualitatively describes the physics of these 
systems and provides a useful language to discuss their 
behavior. However, it quantitatively disagrees with the 
Boltzmann equation to which it is an approximation. 

The differences between the results found from the 
Boltzmann equation and those from the drift-diffusion 
calculation arise for a couple of reasons. One reason is 
related to the failure of the drift-diffusion calculation in 
other cases of in-plane transport. While currents consist 
of electrons moving in all directions, in the drift-diffusion 
approximation, that motion in all (three-dimensional) di- 
rections is averaged, leaving a single direction for the cur- 
rent. A consequence of this averaging is that the model 
misses the injection of current and spin current moving 
parallel to the interface from one layer to the other. In 
the case of CIP GMR, the lack of injected spin currents 
flowing from one layer to the other eliminates any CIP 
GMR. In the present case, the approximation misses the 
injection of parallel current across the interfaces. The in- 
jection (or reduction) of the current flow in the plane of 
the interface can change the resulting torque by a factor 
of two or more. Another source of disagreement between 



the calculations is that the distribution function has a 
different form than that assumed when formulating the 
boundary conditions in the magnetoelectronic circuit the- 
ory. The differences between the assumed distribution 
function and that near the interface in the Boltzmann 
equation gives quantitative differences between the two. 

These models provide a framework that naturally in- 
cludes both the torques due to the bulk spin Hall ef- 
fect and the spin transfer torque and the torques due to 
the interfacial spin-orbit coupling. These two torques 
are the current induced torques that arise from spin- 
orbit coupling and which are independent of the gra- 
dient of the magnetization. The models we treat are 
qualitatively similar to previous models for the spin-Hall- 
induced torque but differ substantially from the models 
used to describe the Rashba torque. Those latter models 
are based on a two dimensional treatment of the trans- 
port that gives rise to a current induced spin accumula- 
tion and a predominantly field-like torque. The Boltz- 
mann equation approach includes the interfacial spin- 
orbit coupling in the boundary conditions of a three di- 
mensional transport calculation. Nevertheless, this ap- 
proach gives very similar results to the two-dimensional 
calculations. The interfacial spin-orbit coupling gives rise 
to a torque that is predominantly field-like. Depending 
on the specific parameters appropriate for a particular 
system, either the field-like torque or the damping-like 
torque may dominate. 

While the experimental situation is still controversial, 
there is experimental evidence that both a damping-like 
torque and in some systems a field-like torque play an im- 
portant role in the dynamics. The model developed here 
through the Boltzmann equation captures the physics for 
both. Unfortunately, it is difficult to make the model 
predictive rather than explanatory. While many of the 
transport parameters are known for thick films, they are 
likely to change significantly in thin films. In fact, many 
vary with varying thicknesses of the films. This model 
does capture the variation with the thickness of the non- 
magnetic layer, but does not describe the rapid vari- 
ation with ferromagnetic film thickness found in some 
systems. This behavior, coupled with the variation of 
behavior with the order of growth, suggests to us that 
structural aspects of the samples vary with thickness or 
growth order. Examples of process that might contribute 
to this variation include, intcrdiffusion, strain relief, or 
grain size. 
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Appendix A: Boltzmann Equation 

The Boltzmann equation is a semiclassical approach 
based on the approximation that in some small but not 
too small region of space it is possible to define electron 
wave packets that have both a well defined momentum 
and a well defined position. It is related to a density ma- 
trix approach which neglects all of the coherence between 
states with different wave vectors. The basic quantity of 
interest is the distribution function, /(k, r), which is the 
probability to find an electron with wave vector k at po- 
sition r. 

The straightforward generalization of the Boltzmann 
equation for spin polarized systems is to have separate 
distribution functions for up and down electrons, f' and 
/■k This is the approach used by Camley and Barnas 60 
to model GMR. For systems with spin-orbit coupling, in 
which spins can point in arbitrary directions, there are 
two related approaches to generalizing the distribution 
function. In analogy with the density matrix, the distri- 
bution function can be generalized to a 2x2 Hermitian 
matrix in spin space, f. Alternatively, the same informa- 
tion can be captured by four real distribution functions 
f a , a — 0, x, y, z related to f by 



f a = Tr[<7 a f], 



where 



(To 



1 
1 

1 

1 

-i 

1 



(Al) 

(A2) 
(A3) 
(A4) 



1 

-1 

{(T x ,(Ty, (T z ) 



(A5) 
(A6) 



The generalized Boltzmann equation is 



dfg drdfa 
dt dt dr 

— ^ [fp(t,T,k,n)], 
at coil 



dk dfg 
dt dk 



- jH^fjCap-f 



(A7) 



where the collision term on the right hand side depends 
on all four distribution functions. The last term on the 
left hand side describes spin precession in ferromagnetic 
layers, where the electron spins precess in the exchange 
field H cx . The time derivatives of r and k are given by 



dr 

~dt 

dk 
i— - 

dt 



V k ,n 

-eE, 



(A8) 
(A9) 



where vt, n is the velocity of the electron and E is the 
electric field. For the linearized Boltzmann equation, we 
make the replacement 

/„(k) -► / eq (e(k))5 a , + ffa (K)/^(e(k)), (A10) 

Upper case K refers to wave vectors restricted to the 
Fermi surface. f^ q is the energy derivative of the Fermi 
function, emulating a Taylor series expansion of the dis- 
tribution function around equilibrium. 



After some standard algebra, the linearized Boltzmann equation can be cast into the form 

= -R a , a ,(K i )g a ,(K i )+ f dK fP a . a > (K< , Kf)g a > (K/) (All) 

Jfs 

The first term on the right hand side is the scattering out term and the second term is the scattering in term. The 
former describes collision processes that reduce the occupancy of a state and the latter those that increase it. 



VK 



dg a (Ki) 
dr 



eE • v K <5 Q ,o + ^H'fg 1 {K)e a p 1 



In the ferromagnet, where we neglect spin-orbit coupling, the scattering is diagonal in a coordinate system aligned 
with the magnetization. For the magnetization pointing in a general direction (sin# cos </>, sin# sin0, cos 9), the spin- 
dependent scattering matrix for the scattering out terms is 

R a , a ' = U^pRp^Sp^'Up/rf. (A12) 

The diagonal scattering matrix is 

i? diag = (R\R,R,R l ), (A13) 

in terms of the majority and minority scattering rate RJ = \/t* and R^ — and the transverse scattering rate 
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is taken to be the geometric mean of the spin-dependent scattering R = 1/ \J t^t^. The transformation matrix is 



U 



( 1/V2 



\ 1/V2 




(A14) 



The first matrix transforms between the majority/minority description on one hand and the Cartesian description 
that is used in the rest of the calculations. The second two matrices rotate the coordinate system. 



The scattering in terms are similar. For these terms, 
where spin-orbit coupling does not play a role, P a , a ' is 
independent of wave vector and equal to R a , a '- There is 
an additional contribution from spin flip scattering of the 
form R S a , a > (1 — <5 Q ,o) in terms of the spin flip scattering 
rate R sf — l/r s f. The last factor restricts the scatter- 
ing to the spin distribution functions and not the charge 
function. 

In the non-magnet, the spin independent scattering is 
included through a term of the form i? N 5 ajQ ,/ in terms 
of the scattering rate i? N = 1/r, and spin-flip scattering 
through a term of the form i? Ns <5 QiQ '(l — <5 Q ,o) m terms 
of the spin flip scattering rate i? Nsf = l/r s f. 

Spin-orbit scattering is more complicated than the 
scattering processes described above because the scatter- 
ing rates depends on the initial and final momenta. Engel 
et alM- give the contribution to the collision integral as 



rff 

dt sH 



nihk F ^ da 



77T " ' dfl 

k f 



Ehk 
— 
m* 

k f 



(A15) 



■mS(<p)a ■ Ij-^t/oCkO + /o(k/)] 



where ip is the angle between and k^. This form is 
based on assuming that the spin-orbit scattering is weak 
and only keeping quantities lowest order in the spin-flip 
scattering. We follow this approximation but also in- 
clude the scattering that gives the inverse spin Hall ef- 
fect. The complete story story is much more complicated 
because there is also scattering between different spin 
channels. We assume that I(ip) is a constant and that 
S(ip) = S\ki x k/|, where S is now a constant. Then, in 
our notation, we have 



df a 
dt sH 



hk 



E-^( k *)-^ k /)] 



(A16) 



- ISn a [fo(ki) + /o(k/)] 

+ IS8 afi n a > [f a > (kj) + f a * (k f )} 

^[-/5n Q /o(k/) + IS5 afl n a ,f a ,(]i.f)} 



taming S a fi are the additional terms that give the inverse 
spin Hall effect. In the second step, we have dropped the 
isotropic part because it is simply another contribution 
to the isotropic scattering. We also find that for the 
remaining parts, the scattering out contribution is zero 
because Ylf n a — 0. Finally, we absorbed the velocity 
and the impurity density factor into the scattering rate. 
Translating to the notation we have been using gives 



Pn 



PS[k i pkf 1 

} (A17) 



where n = k, x k 



kjfikt-t. The terms con- 



Given the scattering matrices, we find the general so- 
lutions of the Boltzmann equation in each layer using 
the techniques described in Ref. [fH These are matched 
together at the interface through boundary conditions 
based on Eq. (flUl) and Eq. (TT21) and subjected to diffuse 
or specular boundary conditions at the out interfaces. 

The default values of the parameters we use are given 
in Table HI These have been chosen to approximately have 
values appropriate for Co/Pt bilayers with vacuum on ei- 
ther side. These parameters are either input parameters, 
calculated numerically based on the input parameters or 
determined analytically from them. In the last case, the 
evaluated expression is given in the table. 



Appendix B: Analytical expression based on circuit 
theory 

In this Section, we present an analytical expression 
of spin torque caused by the spin Hall effect (only) in 
NM|FM bilayer structures where NM has strong spin- 
orbit coupling and is thus subject to the spin Hall effect. 
The derivation is based on the drift-diffusion model 64 
with the boundary condition of the magnetoelectronic 
circuit theory*^ This result does not account for Rashba- 
like interactions at the interface but provides a closed- 
form result in their absence. 

The bulk transport is determined by Eqs.([2]|8]) and 
the current across the interface is calculated using mag- 
netoelectronic circuit theory i 67 i 73 The charge and spin 
currents satisfy the boundary conditions at the NM|FM 
boundary (z — 0) given by 

j z = (G t + G^A/i - (G t - G l )Afj, s ■ M 
z Q = Re[G n ](2A/x s x M x M) - Im[G n ](2A/x s x M) 
-(G t + G i )A/x s + (G t -G^)A/iM, (Bl) 
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TABLE I: Default parameter values. Parameters for the fer- 
romagnet (F), are chosen to be roughly those for Co as in 
Ref. UA an d those for the non-magnet (N) , to be roughly those 
for Pt as in Ref. [3^ . A = dpt is a mean free path and £ si is a 
spin diffusion length. The rest of the parameters are defined 



in the text . 








An 


input 


2.43 


imi 


\ 8 f 
/X x 


input 


1 A 7 


nm 


^sH 


input 


ll.o 


imi 


\t 
N 


input 


1 P. OK 


nm 


\i 
N 


input 


6.01 


nm 


\sf 
/X F 


input 


oocn 


nm 


\ 


input 


U.zoo 


nm 


fcp 


input 


16 


nnf 1 


ItQ 


input 


0.42645 






input 


0.20055 






input 


0.04 




Osa 


computed 


-0.059 




esf 
f-N 


computed 


2.57 


nm 


«F 


computed 


69.3 


nm 


TN 


e 2 2A N . 2 
h 3^ nK F 


0.005 


nm^Q -1 


Pn 


1/(T N 

e 2 A^+A^ , 2 
ft 3tt 2 "^F 


20 


/ificm 


<rp 


0.02 


nm _1 ri _1 


Pf 


1/ctf 


5 


/if2cm 


G , t 


evaluated 


6.66xl0 14 




G l 


evaluated 


3.96xl0 14 




Re[G n ] 


evaluated 


5.94xl0 14 


n-hn- 2 


Im[G^] 


evaluated 


0.86xl0 14 


n-hn- 2 



where G^ and G^ are interface conductances for majority 
and minority spins, aligned antiparallel and parallel to M 
respectively, G'^ is the mixing conductance, A/x = fi(z = 
+0) — [i(z = —0) is the chemical potential drop over the 
interface, and A/x s is the spin chemical potential drop 
across the interface. 



By solving the bulk equations, Eqs.([5][S]), the boundary 
conditions at the interface between the materials, Eq. 
(|B1[) . with the additional boundary conditions, j z = at 
z = 0, and j s = at z = +tp and z — —t^, where tp 
and <n are the thicknesses of FM and NM, respectively, 
one obtains coefficients for two vector components of spin 
torque as in Eq. ([T3]) 



^ 2f „ 2Im[G^] 2 (l + 4)ff + ^[^J [2(l + 4)i>[^] + (l-«M 
[2Im[Gft4-](l + 4 )^] 2 + [2(1 + 4)Z£Re[Gt+] + (i _ 4 )(JN ] 2 

r f = 0SH 2( KN - l) 2 ^ M^](l-^ )g K (B3) 

[2Im[Gti](l + « N )^] 2 + [2(1 + 4)Z£Re[GT+] + (1 - ^ )a N ] 2 

where #sh = ctsh/cn is the spin Hall angle, on is the conductivity of NM, kn = exp(— t-^/l^), and iif is the spin 
diffusion length of NM. We note that both Td and T{ do not depend on the magnetization direction for the case with 
the spin Hall effect only. Thus for the case with the spin Hall effect only, the angular dependence of the torque on 
the magnetization direction is completely determined by the cross products in Eq. (| 14[) . These damping and field like 
torques are plotted in Fig. (j6|) (dashed lines). These results are independent of the thickness of the ferromagnetic 
layer because the boundary conditions force the transverse spin current to be zero in the ferromagnet. All of the 
dependence of the non-magnetic layer thickness is captured in the factors of kn • 
The mixing conductance in Eq. (|B2[) and Eq. (IB3I) is evaluated by the integral 

'■"--'jlj^' 1 -'-'- 1 m 

FS 

Here, "up" and "down" spins are defined with respect to the magnetization, so that s = 1 corresponds to a minority 
electron with its spin parallel and moment antiparallel to the magnetization. FS refers to integrating over the Fermi 
surface. The mixing conductance becomes 

G,i = x/(0< 1 -^'- (B5 » 

FS 
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On the right hand side, f & n d i refer to majority and minority electrons respectively. For the model treated here, 
with a delta function sheet potential, Eq. (j!2[) gives the reflection amplitudes, so the mixing conductance becomes 



2nh 



dxx 1 



ix — —ix — 



1 



2nh 1 2 2(ut + u 4-) 



2(ttt + u+) 



itHn 



1 + u± 



^ln 



1 + u f 



(tt - 2tan~ 1 (u ; )) - u f (7r - 2tan~ 1 (u T ))] 



(B6) 



Here = uq — u cx and = uq + u cx . If > u} as is the case here, both the real and the imaginary parts of the 
mixing conductance are positive. 
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